Detecting a Majorana-Fermion Zero Mode Using a Quantum Dot 
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We propose an experimental setup for detecting a Majorana zero mode consisting of a spinless quantum 
dot coupled to the end of a p-wave superconducting nanowire. The Majorana bound state at the end of the 
wire strongly influences the conductance through the quantum dot: driving the wire through the topological 
phase transition causes a sharp jump in the conductance by a factor of 1/2. In the topological phase, the zero 
temperature peak value of the dot conductance (i.e. when the dot is on resonance and symmetrically coupled to 
the leads) is e^/2/i. In contrast, if the wire is in its trivial phase, the conductance peak value is e^/h, or if a 
regular fermionic zero mode occurs on the end of the wire, the conductance is 0. The system can also be used 
to tune Flensberg's qubit system [PRL 106, 090503 (2011)] to the required degeneracy point. 

PACS numbers: 73.21. -b, 74.78.Na, 73.63.-b, 03.67.Lx 
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Majorana fermions, an exotic type of quasi-particle with 
non-Abelian statistics, are attracting a great deal of attention 
due to both their fundamental interest and their potential ap- 
plication for decoherence-free quantum computation. Several 
ways to realize unpaired Majorana fermions in a vortex core 
in a /?-wave superconductor^"^ and superfluid^'^ have been 
proposed. Majorana bound states (MBS) may also be real- 
ized at the ends of a one-dimensional;?- wave superconductor- 
for which the proposed system is a semiconductor nanowire 
with Rashba spin-orbit interaction to which both a magnetic 
field and proximity-induced 5-- wave pairing are added. ^^'^^ In 
view of these proposals, how to detect and verify the exis- 
tence of MBS becomes a key issue. Suggestions include noise 
measurements r^ii^ resonant Andreev reflection by an STM^ 
and 47r periodic Majorana- Josephson currents i^'iJ^ 

With regard to quantum computation, the braiding of Ma- 
jorana bound states in a network of wires by applying a 
"keyboard" of individually tunable gatesi^ leads to non-trivial 
computation. Note that all the detecting methods proposed 
to date,^~^^ involving electron transfer into or out of MBS, 
will destroy the qubit information. In addition, such braid- 
ing can not result in universal quantum computation; it must 
be supplemented by a topologically unprotected 7r/8 phase 
gate J^ Recently, Flensberg introduced a system consisting of 
a quantum dot coupled to two MBS (MBS-dot-MBS) through 
which this 7r/8 phase gate can be achieved. ^^ A key point is 
that the system must be fine tuned so that the ground state is 
degenerated^ 

In this work, we consider a spinless quantum dot coupled to 
a MBS at the end of a/?- wave superconducting (SC) nanowire, 
and study the conductance, G, through the dot by adding two 
external leads (schematic in Fig.[T]). We find that the con- 
ductance is independent of the properties of the MBS, the 
nanowire, or the superconductor. The dependence of G on 
the dot properties has the same functional form whether an 
MBS is present or not. Therefore, the conductance behavior 
can be conveniently summarized by its peak value, when the 
dot is on resonance and symmetrically coupled to the probing 
leads. It is e^ /2h in the topological SC phase, Gpeak = 1/2, 
in contrast to that for a dot coupled to a regular fermionic zero 
mode, Gpeak = 0, as well as to that for a dot coupled to the 
wire in its topologically trivial phase, Gpeak = 1- Thus, as 
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FIG. 1: (color online) (a) Sketch of dot-MBS syste:the semiconduc- 
tor wire on a 5-wave superconductor surface, and a magnetic field 
perpendicular to the surface {z direction). The dot couples to one 
end of the wire; the conductance through the dot is measured by 
adding two external leads, (b) Majorana chain representation for 
leads-dot-MBS system (Cpeak = e^/2/i). (c) Dot-leads system with 
nothing side-coupled (left) and Majorana chain representation (right) 
(Cpeak = e^ /h). (d) Dot-leads system with side-coupled regular 
fermionic zero mode (left) and Majorana chain representation (right) 

(G'peak = 0). 

the wire is driven through the topological phase transition, the 
conductance shows a sharp jump by a factor of 1/2. The con- 
ductance through the dot is, then, a probe of the presence of 
the MBS. Note that direct transfer between the MBS and dot 
is not necessary, though dephasing of the qubit is introduced 
when the dot is on-resonance. Such a "less invasive" sens- 
ing method provides a potential way to probe a MBS without 
totally destroying the information in the qubit. We also con- 
sider coupling the dot to both ends of the wire (two MBS), 
with a magnetic flux $ through the loop. The conductance as 
a function of phase shows peaks at <l> = (2n + l)7r<l>o which 
can be used to tune Flensberg's qubit systen>iS to the energy 
degeneracy point. 

Single MBS — We consider the setup shown in Fig. (Ha) in 
which a spinless quantum dot is coupled to the end of a semi- 
conductor nanowire with strong Rashba spin-orbit interaction, 
proximity-induced ^-wave superconductivity, and a magnetic 
field B}^'^^ We assume the nanowire and superconductor are 
not grounded and have a negligible charging energy. The mag- 
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The standard equation of motion method yields an exact ex- 
pression for the Green function^^ 



FIG. 2: (color online) Spectral function of the quantum dot in the 
on-resonance (cd = 0) and symmetric (Tl = ^r — r/2) case, (a) 
Coupling from dot to MBS (A) and leads (F) varies at fixed tu — 0. 
Solid lines: F = 0.2 and A from to 0.1. Dashed lines: A = 0.02 
and F from 0.05 to 0.1. The spectral function evolves from a simple 
resonant tunneling form in the absence of coupling to a three-peak 
structure; the middle peak is a direct result of the Majorana zero 
mode, (b) MBS-MBS coupling strength varies at fixed F = 0.2, 
A = 0.1. Note that A(uj — 0) = 1/2 whenever a Majorana is 
coupled. The unit is chosen so that the lead band width is Dl = 40 
for all calculations. 

netic field is smaller than the superconductor's upper criti- 
cal field, but the Zeeman splitting Vz = gfiBB/2 must be 
large enough for the wire to be in the topological SC phase, 
Vz > yA^+T^ where A is the SC order parameter and fi is 
the chemical potential of the wire. Isolated Majorana fermion 
zero modes r]i and r]2 appear in this case at the two ends of 
the wire. Suppose the dot is coupled to r]i and the operators 
d^ (cl^) create an electron in the dot (leads). The Hamiltonian 
can then be written a^ 

H = i^Leads + ^Dot + ^D-L + i^MVlV2 + A(d - d'^)r]ii W 

where il^Leads = E/c T.a^L,R ^kc\^Cka describes the left and 
right metallic leads with chemical potential /iiead = 0, i^oot = 
Cdd'^d describes the dot with a gate tunable level Cd, and 
i^D-L = Ea=L,i? Y^k '^a(4a^ + ^'^') ^escribes the coupling 
between the dot and the leads, cm ^ e"^/*^ is the coupling 
between the two Majorana bound states, where L is the length 
of the wire and ^ is the superconducting coherence length. 

The last part of H describes the coupling between the dot 
and MBS. Here, we assume that the Zeeman splitting is the 
largest scale Vz ^ I'^biasl, ^, T, A, where Vbias is the source- 
drain voltage, T is temperature, and F = F^ + F^^ is the 
dot-leads coupling with To, = TT\Va\^ pQ and po the density 
of states of the leads. In this case, one need only consider a 
spinless single level in the dot. It is helpful to switch from the 
Majorana fermion representation to the completely equivalent 
regular fermion one by defining rfi = (/ + /^) /\/2 and 7/2 = 
i {f — f^) /V^- The last two terms in H become 

Hubs = eM{Pf - ^) + A(d - d^) (/ + p) /V2 . (2) 

The linear conductance through the lead/dot/lead system is 
related to the Green function of the dot level, Gf^(cj), by 
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For Cm = and e^ = 0, one has G^^{uj ^ 0) = l/2{uj + 
zF), and so the on-resonance (e^ = 0) and symmetric (Vl = 
Vr),!.^. peak, conductance at zero temperature is 

Gpeak = -{eyh)rim[G^a{u ^ 0)] = ey2h . (6) 

This result is distinct from both the case of a dot coupled to a 
regular fermionic zero mode, which gives Gpeak = 0,^^ and that 
of a dot disconnected from the wire, for which Gpeak = e'^/h. 
For asymmetric coupling (Vl 7^ Vr), there is a pre-f actor 
"^^L^n/i^L + Ti?)^ for all cases. Therefore, the signature 
of the Majorana fermion is that the conductance is reduced by 
a factor of 1/2. 

To further understand this result, we rewrite the model in 
the Majorana representation.^ The probe leads are described 
by two semi-infinite tight-binding fermionic chains Ci (i = 
..., — 1, 0, 1, 2, ...) joined at the dot, i = 0. By transforming 
to the Majoranas (Greek letters) , (3i = {ci + cl)/\^ and 
7i = {—ici + icl)/\^, our model reduces to two decoupled 
Majorana chains, as shown in Fig.[Hb). The side-coupled 
MBS in the lower chain corresponds to the MBS rji . The con- 
ductance through the dot is, then, the sum of the conductance 
from two decoupled Majorana chain G = G^pp^' + G^^^^^ 

Consider now two other cases. First, for a system with- 
out a side-coupled mode, the Majorana representation leads 
to two decoupled chains as shown in Fig.IHc). Second, for 
a system with a side-coupled regular fermionic zero mode, 
the Majorana representation consists of two decoupled chains, 
each of which has a side-coupled MBS [Fig.fHd)]. For both 

cases, ilf^PP^^ = -ilflower^ ^^^ ^j^^g g^upper _ flower 3^^^^^ 

the peak conductance for a dot with (without) a side-coupled 
regular fermionic zero mode is (e^//i), the result for a sin- 
gle Majorana chain with (without) a side-coupled MBS is 
(e^/2/i).Therefore, the conductance of our model [Fig.[TJb)] 
is Gpeak = + e^/2h = e^/2h. 

The spectral function ofthe dot, A (cj) = -2FIm[Gf^(a;)], 
is shown in Fig.[2ta) for several values of the dot-MBS cou- 
pling A and dot-lead coupling F for eM = 0. The energy unit 
is chosen so that the lead band width is Dl = 40 through- 
out the paper. Consistent with our assumption that the Zee- 
man splitting is the largest energy scale, we consider the spec- 
trum for only the spin-down channel. For A = 0, the spec- 
tral function reduces to the result of the resonant level model. 
For small dot-MBS coupling (A = 0.02,0.05), the spectrum 
shows two peaks at a; ^ ±A which come from the energy 
level splitting caused by coupling to the MBS. As we increase 
A with fixed F = 0.2, the two peak structure evolves into a 
spectrum with three peaks, showing clearly the presence of 



the Majorana zero mode. Note that the zero frequency spec- 
tral function always gives A(cj = 0) = l/2as long as eM = 
and A 7^ 0. For small dot-MBS coupling (A = 0.02), the three 
peak spectrum also appears upon decreasing F. 

The dot spectrum for different strengths of MBS-MBS cou- 
pling Cm appears in Fig.[2tb). Even for very small coupling 
Cm = 0.02, the zero frequency spectrum shows A((jJ = 0) = 
1 not 1/2. The width of the narrow peak is proportional to 
Cm- For large coupling {cm = 0.3), the spectrum reduces 
to the resonant level result along with two additional small 
peaks at a; ~ ±eM corresponding to the energy of the effec- 
tive Dirac fermionic state /. If the wire is long enough so that 
Cm <C T, a, one can still observe the Gpeak = e^ /2h signature. 

More Realistic Wire — To analyze the robustness of the 
MBS signature in the real physical system, the single MBS 
in Eq. ([T]) is replaced by the whole nanowire^^'^^ shown 
in Fig. (Ha). We study numerically a lattice tight-binding 
Hamiltonian,^^ il^wire = i^o + i^Rashba + ^sc, where Hq in- 
cludes nearest-neighbor hopping along the wire {y direction), 
a chemical potential leading to half filling (/i = 0), and a mag- 
netic field perpendicular to the surface (i direction) causing 
the Zeeman splitting Vz . The Rashba spin-orbit interaction is 
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-i aRw\^^ z ■ (a^,/ X x) w- ' + h.c. (7) 



where wl ^ creates an electron with spin index s on site i of the 
wire and a are the Pauli matrices. Finally, the ^-wave pairing 
term with superconducting order parameter A is 



i7sc=A^ 



wl^w, 



t 



h.c. 



(8) 



The Bogoliubov-deGennes equation is constructed from i^wire 
by the standard Nambu spinor representation (including the 
same Zeeman splitting Vz in the dot) and then solved by a 
recursive Green function method.^^'^^ 

The dot spectral function is shown in Fig. [3] for several 
values of the SC order parameter A and Rashba interac- 
tion strength aR (for an on-resonance, symmetrically cou- 
pled dot)i21 When the wire is in the topologically trivial phase 
(A > Vz, no MBS), the spectrum is similar to the resonant 
level result [Fig.Oa)]. In contrast, when the wire is in the 
topological SC phase (A < 1^ = 6, /i = 0), the value of 
the spectral function at zero frequency is 1/2. For A small 
(A = 0.5), the spectrum shows two peaks, but upon increasing 
A (A = 1,3), the two peaks become more separate and the 
three-peak structure emerges. Curiously, a further increase 
of A (A = 4.5, 5.2) leads to a smaller separation between 
the outer peaks. Similar phenomena occur upon varying the 
Rashba interaction aR. increasing aR leads to first an in- 
crease in the splitting of the outer peaks {aR = 1, 4, 10) and 
then a decrease {aR = 15, 25). 

The non-monotonic shifts in the positions of the outer peaks 
can be understood as follows. When A or aR is small, the/?- 
wave SC pairing fp is weak, leading to a less robust MBS and 
small peak splitting. On the other hand, when A is large and 
close to the transition value A = 14, SC pairing between the 
lower and upper band^ makes the MBS less robust. For large 
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FIG. 3: (color online) Dot spectral function and peak conductance in 
the more realistic nanowire case (the dot is on-resonance and sym- 
metrically coupled to the probe leads). A{uj) for different values of 
(a) the SC order parameter at fixed aR = 2 and (b) the Rashba inter- 
action strength at fixed A = 3. The results are qualitatively similar 
to those of the simple model (Fig. 2). (Parameters: F = 0.1, A = 0.3, 
and Vz = Q.) (c),(d) Conductance as a function of Zeeman energy for 
different temperatures at fixed A = 3. The sharp change at 14 = A 
is a signature of the topological phase transition. (Parameters: (c) 
aR^2, A = 0.1, r = 0.1; (d) ^^ = 10, A = 0.3, r = 0.08.) Through- 
out, /i = 0, the hopping in the nanowire t = 10 corresponds to a band 
width D — 40, and the wire consists of 1000 sites. 



aR, the eigenfunction of the lower band at the fermi surface 
has a large spin-up component, while the dot and leads are 
spin-down due to the Zeeman splitting; therefore, the coupling 
between the dot and MBS is suppressed. As a function of both 
parameters, then, there is non-monotonic behavior. 

To detect the MBS, a clear signature appears in the con- 
ductance as a function of Zeeman splitting [Fig.^c) and (d)]: 
the conductance at zero temperature shows a sharp jump at 
14 = A due to a topological phase transition.— For 14 < A, 
the wire is in the topologically trivial phase, and the peak con- 
ductance is c^ /h. For Vz > A, the wire is in the topological 
SC phase in which a MBS appears, and the peak conductance 
is c^ /2h. (Both of these values are multiplied by the factor 
^^L^R/iXh + Ti^)^ for asymmetric coupling to the leads.) 
At finite temperature, the jump becomes a crossover, which 
is still quite sharp near the transition point. For small aR, A 
and large F, the spectrum has two peaks, so the finite T con- 
ductance is larger than c^ /2h [Fig.^c)]. For large aR, A and 
small r, the spectrum has three peaks, causing the finite T 
conductance to be smaller than c^ /2h [Fig.Od)]. 

We emphasize that the change in conductance by a factor of 
1/2 is universal as long as the MBS appears and couples to the 
dot. With regard to the effect of disorder in the wire,^^ a short 
range impurity potential does not affect the MBS and thus the 
G'peak = c^ /2h result, while a sufficiently strong long range 
impurity potential may induce mixing of the MBS at the two 
ends and therefore lead to Gpeak = c^ /h as shown in Fig.[2b). 

Two MBS — Consider the geometry proposed by 
Flensberg^^ for implementing a 7r/8 phase gate: a dot 
coupled to both ends of the nanowire — and hence to two 
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FIG. 4: (color online) (a) Conductance for MBS-dot-MBS system 
as a function of the phase = $/$o for different temperatures; 
the dot is on-resonance and symmetrically coupled to the STM tips. 
(T = 0.01, 0.005, 0.00125, and 0, from top to bottom; parameters 
are Ai = A2 = Fi = r2 = 0.1.) This curve does not depend 
on the value of IA1/A2I. (b) Sketch of MBS-dot-MBS system. The 
two MBS appear at the ends of the nanowire; $ is the magnetic flux 
through the loop. The conductance is measured using dual-tip STM, 
allowing tuning to the degeneracy point. 



MBS — with magnetic flux $ through the loop, as shown in 
Fig-Sit)). The conductance through the dot is measured using 
two external leads; since electron tunneling between the dot 
and environment should be avoided during qubit operation, a 
dual-tip STM setup2&22 is proposed so that one can remove 
the external leads after tuning the system. The Hamiltonian 
of this MBS-dot-MBS system^^ can be written as 



H = Cdd^d + {\ld) - \id)r]i + i{\ld^ + \2d)r]2 



(9) 



The phase difference between the two couplings, (j) = 
2 arg(Ai/A2), is related to the flux <l> via (/) = ^/^q, where 
^0 = h/2e. Without loss of generality, we take Ai to be real 
(Ai = I All and A2 = |A2|e~*'^/^), and the Hamiltonian re- 
duces to H = ed^d-\- X{d'^r]i2 + ^i2<^) where r]i2 = {\Xi\r]i + 



ze^^/2|A2|7?2)/A andA=V|AiP + |A2P. For (/>= (2n + 1)^ 
(n integer), we have 7/12 = 7?|2. In this case, the dot is ef- 
fectively coupled to a single MBS 7^12; therefore, the T = 
on-resonance conductance is e^ /2h. For (j) ^ (2n + l)7r, 
we have 7/12 7^ r/12 corresponding to a regular fermionic zero 
mode, for which the T = on-resonance conductance is zero. 
Following the method for single MBS, one can exactly 
solve for the dot Green function G^^{uj) in this two MBS 
problem in the case of eM =0: 

GS/H = {[Gf°H]-' - ^(c.) - B(c.)}"' (10) 
where A{uj) = -iT + (|Ai \^ + |A2p)/2w and 



B{w) = 



1 

4u;2 



[|Ai|4 + |A2|^ + 2|Aip|A2pcos((/>-^)] 



c^ + erf + ir-(|Ai|2 + |A2P)/2c^ 



(11) 

The conductance peak value as a function of the phase 
difference ^ can be obtained from Eq. ([3]) and is shown in 
Fig.lUa). For T = 0, the G = e^/2/i signature appears only 
at (/) = (2n + l)7r, corresponding to the energetically degen- 
erate state in Flensberg's qubit, ^^ with G = otherwise. For 
T 7^ 0, the peak width becomes finite; note that the peak is 
fairly wide even for T = 0.01 but that the temperature is still 
low enough to see the MBS. By tuning the conductance to a 
resonance peak, one can tune the MBS-dot-MBS to the de- 
sired degenerate energy point. 
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